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43 ■ 1 Introduction 

> '■ 

The discovery of the Higgs boson is one of the main goals of searches at the present 
and next hadronic colliders, the Tevatron and the LHC. 

One appealing production mode, the double pomeron exchange (DPE) one, was 
proposed some time ago in Refs. [1,2]. In the following papers this subject was 
discussed from different perspectives [3-11]. Despite some progress the serious un- 
certainties are still present that do not allow to get fully reliable predictions needed 
for future experiments. This reflects our present limitted understanding of the na- 
ture of the diffractive (pomeron) reactions. 

The best way to reduce these uncertainties is to study other double pomeron ex- 
change processes and compare them with existing data. A particularly enlightening 
process is the DPE production of two jets (dijets). Such a process was originally 
discussed at the Born level in [12]. Later the dijet production was studied in [5, 13] 
and in [8-11,14-19]. 
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One generally considers two types of DPE events when colliding hadrons remain 
intact, namely exclusive and central inclusive one (or central inelastic). In the 
exclusive DPE event the central object H is produced alone, separated from the 
outgoing hadrons by rapidity gaps: 

pp — > p + gap + H + gap +p. (1) 

In the central inclusive DPE event there is an additional radiation X accompanying 
the central object H: 

pp — > p + gap + HX + gap + p. (2) 

Recently, using the Bialas-Landshoff [2] model for central inclusive double diffrac- 
tive production the cross-section for gluon jet production was calculated [18, 19]. In 
this model in some approximation pomeron exchange corresponds to the exchange 
of a pair of non-perturbative gluons which takes place between a pair of colliding 
quarks [20]. The obtained results together with those for quark- ant iquark jets cal- 
culated some time ago [15] give the full cross-section for dijet production in double 
pomeron exchange reactions. The model was found [19] to give correct order of 
magnitude for the measured [21] central inclusive dijet cross sections. 

In this Letter we propose the simple extension of this model to the exclusive 
processes. We find the extended model to be consistent with the CDF Run I, 
II upper limits [21, 22] for double diffractive exclusive dijet production. We also 
present the predictions for the exclusive Higgs boson production cross sections at 
the Tevatron and the LHC energies. 




Figure 1: Production of Higgs boson H, dijet jj, by double pomeron exchange. The 
colliding hadrons remain intact. 
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2 Central inclusive dijet production 



The matrix element for two gluon jet production in the Bialas-Landshoff model is 
given [18] by the s-channel discontinuity of the diagrams shown in Fig. 2. 





Figure 2: Three diagrams contributing to the amplitude of the process of gluon pair 
production by double pomeron exchange. The dashed lines represent the exchange 
of the non-perturbative gluons. 

The square of the matrix element (averaged and summed over spins and polariza- 
tions) is of the form: 

\M PP \ 2 = 81\M qq \ 2 [F (*0 F (t 2 )} 2 , (3) 
where |.M g9 | 2 is the production amplitude squared for colliding quarks 1 : 

jMj = C f 8l- 2a ^8 2 - 2a ™ exp (2/? ( tl + 1 2 )) R 2 . (4) 

K) K) 

Transverse momenta of the produced gluons are denoted by u 1 and u 2 . The constants 
C and R will be defined later, a (t) — 1 + e + a't is the pomeron Regge trajectory 
with e w 0.08, a' = 0.25 GeV" 2 (h, t 2 are defined in Fig. 1). F (t) = exp(At) is the 
nucleon form-factor with A = 2 GeV -2 . 5±, 5 2 are defined as 5i^ 2 = 1 — ^1,2/^1,2 (^1, 
^2, Pi-, P2, are defined in Fig. 1). The factor exp (2(3 (ti +t 2 )) with (3=1 GeV -2 
[23] takes into account the effect of the momentum transfer dependence of the non- 
perturbative gluon propagator given by (p 2 is the Lorentz square of the momentum 
carried by the non-perturbative gluon): 

D(p 2 ) =D ex V (p 2 /r 2 ). (5) 

The constants C and R are defined as: 

R = 9jdQ 2 Q 2 exp (-3Q 2 ) = 1. (7) 

1 This formula is only valid in the limit of « 1 and for small momentum transfer between 
initial and final quarks. 
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Here G and g are the non-perturbative and perturbative quark gluon couplings 
respectively 2 , r is the range of the non-perturbative gluon propagator (5) and Do 
its magnitude at vanishing momentum transfer. From data on the elastic scattering 
of hadrons one infers D G 2 r = 30 GeV -1 and r = 1 GeV. 

The constant R reflects the structure of the loop integral. Q T is the transverse 
momentum carried by each of the three non-perturbative gluons. R was shown 
explicitly in Eq. (4) for the reason which will become clear in the next section. 

Taking into account (4) we obtain the following result for the differential cross- 
section [19]: 



Here Ce = 81C/(16 (2n) ). E T = \u±\ = \u 2 \ is the transverse energy of one of the 
produced gluons. Ay — y ± — y 2 , y — (yi + 2/2) /2 where y i;2 are the rapidities of the 
produced gluons. For completeness it is necessary to say that the rapidities y± t2 are 
connected with Si, 62 and E T in the following way: 



The result (8) does not take gap survival effect (S 2 ap ) into account i.e. the prob- 
ability of the gaps not to be populated by secondaries produced in the soft rescat- 
tering. From [5,24] we expect that for the Tevatron energy it is about 0.05 — 0.1. 
The factor S 2 ap is not a universal number but it depends on the initial energy and 
the particular final state. Theoretical predictions of the survival factor, S 2 ap , can be 
found in Ref. [25]. 

The main uncertainty in the expression (8) is the value of G 2 /4ir (see (6)). It is 
expected to be [26] about 1 but in fact it should be considered only as an order of 
magnitude estimate. 

Let us now make clear the rather ad hoc nature of many of the assumptions 
inherent in the Bialas-Landshoff approach [2] . The predictions of this model depend 
only weakly on energy (~ s 2e ). This is a consequence of the Regge-like dependence 

2 One should note that the non-perturbative quark gluon coupling G does not depend on the 
scale of the process. 




(8) 



Si y/s = E T exp (yi) + E T exp (y 2 ) , 



8 2 y/s = Ej exp (-yi) + E T exp (-y 2 ) . 



(9) 
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implied by Eq. (4). There are some controversies if such assumption is justified. 
In our calculations we assume the exponential form of the non-perturbative gluon 
propagator (5). As was already stated in [2] there is no reason to believe that the true 
form of D is as simple as this. However, we believe that this is not a serious objection 
to our model. In the Bialas-Landshoff approach the produced object (Higgs, dijet 
etc.) is coupled to the non-perturbative gluons via the perturbative coupling g. It 
is not clear and the question of consistency could be addressed. Finally, let us note 
that estimates in the present Letter are based on the basis of the pure forward 
direction. It was first mentioned in [14] that such approach may lead to incorrect 
results. 

Integrating (8) over the CDF Run I kinematical range [21] for the central inclusive 
production of dijets of E T > 7 GeV we obtain [19] the result to be about 70 nb 
(with G 2 /4n = 1 and no S^ ap ), to be compared with the CDF measurement of 43 
nb (43 ± 26 nb). We thus scale our cross section by a factor of 43/70 ~ 0.6, that is: 



S g 2 a P (v^ = 2TeV) 



(G 2 /4tt)' 

This completes the summary of [18] and [19]. 



= 0.6. (10) 



3 Exclusive dijet production - Sudakov factor 

As was already mentioned the calculation presented in the previous section, based 
on the original Bialas-Landshoff model, is a central inclusive one i.e. the radiation 
is present in the central region of the rapidity. 

In order to describe the exclusive processes one has to forbid this radiation. To 
do it we include the Sudakov survival factor T(Q T ,/i) inside the loop integral (7) 
over Q T . The Sudakov factor T(Q T ,fi) is the survival probability that a gluon with 
transverse momentum Q T remains untouched in the evolution up to the hard scale 
/i = Mgg/2 where M gg is the mass of the produced gluons. The function T(Q T ,fi) 
can be calculated as [5]: 



/ r 2 a s [k\ 
T(Q T)/U ) = exp I -J^ — A_ 



dk 2 . 



k 2 



l-A 



dz 



(11) 



Here A = \k T \ / (fj, + \k T \), P gg (z) and P qg (z) (we take q = u,d,s,u,d,s) are the 
GLAP spitting functions. a s is the strong coupling constant 3 . 

3 In the following we take a s at one loop accuracy i.e. a s (g 2 ) = (47r//?o) (1/ In (g 2 /A 2 )) with 
0o = 9 and A = 200 MeV. 
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Taking into account the leading-order contributions [27] to the GLAP splitting 
functions: 



P gg (z) = 6 [z/(l -z) + (l- z)/z + z(l-z) + S(l-z) (11/2 - n f /3)] , 
P qg (z) = [z 2 + (1 - zf] /2, 



(12) 



we obtain: 



I 



l-A 



zP gg (z) dz = -11/2 + 12A - 9A 2 + 4A 3 - 3A 4 /2 - 6 In A, 



o 




l-A 

P qg (z) dz = 1/3 - A/2 + A 2 /2 - A 3 /3. 



(13) 



Now to describe the exclusive processes we use the formula (8) with R 2 = 1 
replaced by R 2 (/i) where R (/i) is defined as 4 : 



The hard scale fx = M gg /2 can be expressed by E T and Ay in the following way: 



Naturally a question of internal consistency arises. Namely, the Sudakov factor 
uses perturbative gluons whilst in our calculations the Born amplitude (4) uses 
non-perturbative gluons. It is not clear what the non-perturbative gluon is and 
the extension of the original Bialas-Landshoff model to the exclusive processes is 
not straightforward. We hope that taking the Sudakov factor in the loop integral 
into account we obtain an approximate insight into exclusive processes. It should 
be emphasized that at present our calculation is a hybrid of perturbative and non- 
perturbative ideas. 

At the end of this section let us notice that the Sudakov factor (11) does not 
depend on the sum of the dijet rapidities y — (yi + 2/2) /2. This together with the 
observation that y 2 « ((\ + (3)/a') 2 = 144 and AE 2 cosh 2 (^)/s = 6^2 « 1 
leads to the conclusion that the differential cross section for DPE exclusive dijets 
production very weekly depends on the sum of the dijet rapidity y. This feature 
agrees with the observation found in Ref. [13]. Moreover the observed power law 
E- 6 - 5 (with R 2 ~ E- 2 - 2 ) is close to the observation of Ref. [7] (~ E- 7 - 3 ). 

4 Notice that fj, > 1.5 GcV is required so that 9 jf 2 dQ 2 T Q 2 T exp(-3<5^) = 1. 



R(ji) = 9 P dQ 2 Q 2 exp (-3Q 2 ) T (Q T , fj) . 



(14) 



H = £ T cosh(— ). 



(15) 
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4 CDF Run I, II upper limits 



The CDF collaboration has presented results on upper limits on exclusive DPE dijet 
production cross sections. 

At Run I (-^/s = 1.8 TeV) [21] the upper bound for exclusive dijets production 
was measured to be 3.7 nb for the kinematic range of 0.035 < 62 = S p < 0.095 
and jets of E T > 7 GeV confined within —4.2 < y < 2.4 and the gap requirement 
2.4 < ?/g a p < 5.9 on the proton side. 

At Run II (y/s = 1.96 TeV) [22] the upper bound for exclusive dijets of E 1 > 10 
GeV [E T > 25 GeV] was measured to be 970 ± 65(stat) ± 272(syst) pb [34 ± 5(stat) 
± 10(syst) pb]. The kinematics is following 5 : 0.03 < 5 2 = 5 P < 0.1, jets are confined 
within —2.5 < y < 2.5, the gap on the proton side is 3.6 < y gap < 7.5. 

It should be noted that in the above experiments the protons were not detected 
and the DPE events were enhanced by a rapidity gap requirement on the proton 
side 6 . 

Integrating 7 (8) over the appropriate kinematical range we obtain the results 
shown in Table 1. The running coupling constant g 2 /4n, appearing in (6), is evalu- 
ated at 2Ef n i.e. 0.15, 0.14, 0.12 for Ef n = 7, 10,25 GeV respectively. The factor 
S 2 / (G 2 /47r) 2 is taken to be 0.6. 



Transverse 


CDF 


Model 


energy 


upper limits 


SU (G 2 /4vr) 2 = 0.6 


Ej > 7 GeV 


3.7 [nb] 


1 [nb] 


E T > 10 GeV 


970 ± 337 [pb] 


300 [pb] 


E T > 25 GeV 


34 ± 15 [pb] 


3 [pb] 



Table 1: Comparison of the CDF upper limits for DPE exclusive dijet production 
with the results obtained in the presented model. 

As can be seen from Table 1 the obtained results are quite satisfactory and 
encouraging. They are also comparable with those obtained in [28,29] or [17]. 

5 We would like to thank K. Goulianos for a correspondence about this point. 

6 In principle the result (8) should by multiplied by a factor (1 — exp[— 2(A + j3 — 

a' In 8\ ) c ^p^mL) ] ) wnere y^p X is the maximum value of the gap. In the present case, y™p X = 5.9 

and 7.5 for Run I and Run II respectively, this factor is close to 1. 
7 Note an identical final state particle phase space factor ^. 
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5 Exclusive Higgs production 



The matrix element for the Higgs production in the Bialas-Landshoff model is given 
[2] by the s-channel discontinuity of the diagram shown in Fig. 3. The Higgs 
coupling is taken to be through a t-quark loop. 



T 
I 




I 

1 



Figure 3: The diagram contributing to the amplitude of the process of Higgs boson 
production by double pomeron exchange. The Higgs coupling is taken to be through 
a t-quark loop. The dashed lines represent the exchange of the non-perturbative 
gluons. 

The square of the matrix element for colliding (anti)protons has the form [2] 8 : 



\M PP \ 2 = BN 2 s 2 8 2 - 2a{tl) 5 2 - 2a{t2) [F (h) F (t 2 )} 2 

xexp(2(3(t 1 + t 2 ))R 2 . (16) 

The constant B is defined as: 

where Gp is the Fermi coupling constant and a s (Mu) is the perturbative coupling 
evaluated at a scale M#. N is a function of M t jMn. For the Higgs mass Mh < 
2M t ks 350 GeV this function is given by [1]: 



N = 6- 



M? M? ( M 2 \ ( M H \ 2 , x 

-6^ 4-^-1 arcsin-f . (18) 



Mfj M 2 H \ M 2 H J \ 2M t 

It turns out that the structure of the loop integral over Q 7 has exactly the same 
form like that for gluon jets case (7). So to describe the exclusive Higgs production 
we take the result (16) and replace R 2 = 1 by R 2 (/i) with R (//) given by the formula 
(14), where \i = M gg /2 -> M H /2. 

Performing the appropriate calculations we find the differential cross section 
da/dy (not presented in [2]) for DPE exclusive Higgs boson production to be in the 

8 Our result for the matrix element differs from the result of Bialas and Landshoff by a factor 
exp (2(3 (t\ + t-2)). The missing factor 2 pointed out in [3] is also taken into account. 
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form: 

da BN 2 ( s \ 2e R 2 (M H /2) 



(19) 

dy A^a' 2 \M 2 H J ((A + (3) /a' - In [M H /y/s]f - y 2 ' 

Here y, \n(M H / (S^^y/s)) < y < ln(55 nax y/s/M H ), is a rapidity of the produced 
Higgs. In the following we use 5™ ax = #f ax = 5 = 0.1. Since y max = \n(8yfs/M H ) 
and (X + P)/a' = 12 the differential cross section (19) very weekly depends on the 
rapidity of the Higgs. So to get the total cross section for DPE exclusive Higgs 
production it is enough to multiply the above result (19) (at y = 0) by a factor 
f ffi - y min = In (s5 2 /M 2 H ) what leads to the final result: 

BN 2 ( s \ 2t R 2 (M H /2) ( s5 2 



" 4% 3 a' 2 \M 2 H J ((A + (5)/a> - In [M H /^]) 2 \M 2 H 

Now we are ready to give our predictions for DPE exclusive Higgs production 
at the Tevatron and the LHC energies. We also compare our results with those 
obtained in a model developed by Khoze, Martin and Ryskin (KMR model). 

In Table 2 the prediction for the Tevatron energy, y/s = 2 TeV, is shown. The 
mass of the Higgs is taken to be 120 GeV and a s (M H ) is about 0.1. We also include 
the a s virtual correction [5, 30] to the gg — > H vertex factor, so-called if-factor to 
be about 1.5. As was discussed earlier we take 5 = 0.1 and assume S 2 ap / {G 2 /An) 2 
to be 0.6. 



v/s = 2 TeV 


KMR model 


Our model 

S 2 J (G 2 /4n) 2 = 0.6 


<7[fb] 


0.06 


0.005 



Table 2: Our result for DPE exclusive Higgs production cross section for the Teva- 
tron energy. Our prediction is about 10 times smaller than the prediction based on 
the KMR model. 

Before we present the prediction for the LHC energy, y^s = 14 TeV, we have to 
take into account the s dependence of the gap survival factor S 2 ap . Following [25] we 
expect that S 2 ap (y/s = 2 TeV) / S 2 ap (y/s = 14 TeV) ~ 0.4 what allows us to assume 
Sgap/ (G 2 /47r) 2 = 0.25. The obtained result for DPE exclusive Higgs (M H = 120 
GeV) production at the LHC energy is presented in Table 3. 

As can be seen from Table 2 and Table 3 the results for DPE exclusive Higgs pro- 
duction are about one order of magnitude smaller than those obtained in the KMR 
model [5, 6] for the Tevatron energy and about two orders of magnitude smaller for 
the LHC energy. It reflects the possible large uncertainties of the presented ap- 
proach and the general fact that the perturbative QCD predictions, on the contrary 
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y/s = 14 TeV 


KMR model 
Sl p = 0.02 


Our model 

SU (G 2 /4vr) 2 = 0.25 


<r[fb] 


2 


0.015 



Table 3: Our result for DPE exclusive Higgs production cross section for the LHC 
energy. A distinct difference, ~ 10 2 , with the KMR model prediction is observed. 

to the non-perturbative two-gluon-exchange-type models, show a strong increase of 
the cross sections with increasing energy. Hopefully a study of the dijets produc- 
tion as a function of the energy will clearly be able to discriminate between the 
perturbative QCD determinations and the non-perturbative model approaches. 
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